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ABSTRACT 

We present a simple, empirically motivated model that simultaneously predicts the evolution of the 
mean size and the comoving mass density of massive (> 10^^ Mq) early-type galaxies from z = 2 
to the present. First we demonstrate that some size evolution of the population can be expected 
simply due to the continuous emergence of early-type galaxies. The Sloan Digital Sky Survey (SDSS) 
data reveal that in the present-day universe more compact early-type galaxies with a given dynamical 
mass have older stellar populations. This implies that with increasing look-back time, the more 
extended galaxies will be more and more absent from the population. In contrast, at a given stellar 
velocity dispersion, SDSS data show that there is no relation between size and age, which implies that 
the velocity dispersion can be used to estimate the epoch at which galaxies stopped forming stars, 
turning into early- type galaxies. Based on this, we define an empirically motivated, redshift-dependent 
velocity dispersion threshold above which galaxies do not form stars at a significant rate, which we 
associate with the transformation into early-type galaxies. Applying this 'formation' criterion to a 
large sample of nearby early- type galaxies, we predict the redshift evolution in the size distribution 
and the comoving mass density. The resulting evolution in the mean size is roughly half of the 
observed evolution. Then we include a prescription for the merger histories of galaxies between the 
'formation' redshift and the present, based on cosmological simulations of the assembly of dark matter 
halos. Such mergers after the transformation into an early-type galaxy are presumably dissipationless 
('dry'), where the increase in size is expected to be approximately proportional to the increase in 
mass. This model successfully reproduces the observed evolution since z ~ 2 in the mean size and 
in the comoving mass density of early-type galaxies with mass M > 10^^ Mq. We conclude that the 
recently measured, substantial size evolution of early-type galaxies can be explained by the combined 
effect of the continuous emergence of galaxies as early types and their subsequent growth through dry 
merging. 
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1. INTRODUCTION 

Over the past five years it has become clear that 
early-type galaxies do not constitute a passively evolv- 
ing population. Their predominantly old stellar popu- 
lations notwithstanding, substantial evolution occurs up 
until the present day. The main evidence for such evo- 
lution is the increase in comoving mass density by a 
factor of about three from z 1 to the present (e.g . , 
iBell et al.l[2004hl : iBrown eTall [2n07l : [Faber et al.l[2n07h . 
Recently, additional evidence for continuing structural 
evolution of the early-type galaxy population was pro- 
vided by the observation that distant {z > 1.5), quiescent 
galaxies have much smaller sizes than early-type galax- 
ies of the same stellar mass in the present-day universe 
(iDaddi et al.' 2005: Truiillo et al. 2006: Zirm et al."2007'; 
Toft et al.li2007: van D okkum et al. 2008; Cim atti e t al. 
20081 iBuitraio et all l2008l ). This surprising result. 



raising questions about the evolutionary connection of 
these compact galaxies with present -day descendants, 
was recently put on firmer footing bv lvan der Wel et al.l 
(|2008i) : using early- type galaxies with accurate dynam- 
ical masses, thus removing the most important system- 
atic problem that may hamper the z ~ 2 results (the 
absolute uncertainty in the mass estimates), they in- 
ferred substantive size evolution, by a factor of ~ 2, at 
a given inass, b etween z = 1 and present. Recently, also 
iBernardil (|2009f) showed that size evolution continues up 



until the present day, and must be a smooth function of 
redshift. Although the remaining un certainties are not 
negligible (see, e.g.. lHopkins et ani2009a ). it has become 
clear that the observed size evolution cannot be explained 
by systematic uncertainties. 

Because halos that collapsed at earlier epochs were 
denser, and gas-rich, dissipative processes were more 
important for galaxy formation, it has been argued 
that galaxies that formed early are sma ller (e.g., 
[Robertson et~all [2(M iKhochfar fc Sii^ [200l . In this 
scenario, because the number of early-type galaxies in- 
creases with cosmic time, and the more recently formed 
galaxies are larger, there is a direct connection between 
the size and comoving mass density evolution in the 
early-type galaxy population. This is true even without 
individual galaxies changing in size over time. 

This paper explores the extent to which the growth of 
the population can account for the observed size evo- 
lution, and/or whether additional mechanisms to in- 
crease the size evolution of individual galaxies are re- 
quired. As a proof of concept, we first investigate the 
present-day early-type galaxy population for clues that 
indicate that evolution in the mass-size relation may be 
expected purely due to evolution of the population, not 
size growth of individual galaxies. Are small/compact 
early-type galaxies older in terms of their stellar popu- 
lations than large/extended early- type galaxies with the 
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same mass? In Section [3] we show that this is indeed 
the case. Is there a parameter for which age and size 
are independen t? Continuing along the track of previ- 
ous work (e.g.. [Bernardi et al.ll2005t IChang et al.l 120061 : 
iGraves et al. l2008h . which demonstrated that stellar ve- 
locity dispersion (cr*) is the fundamental parameter be- 
hind well-known correlations such as the color-magnitude 
and mass-metallicity relations, we show that a* also cor- 
relates with, and perhaps drives, the age of the stellar 
population. This suggests that formation epoch and cr, 
are related. 

We cast this idea as an empirical relationship between 
the velocity dispersion of an early-type galaxy and the 
redshift zet at which it emerged as an early-type galaxy. 
This model to describe the population evolution of early- 
type galaxies allows us to predict which portion of the 
present-day early-type galaxy population already existed 
at a given redshift z. In turn, such a cr* — zet relation- 
ship in conjunction with the present-day cr* distribution 
thus describes the evolution with redshift of the proper- 
ties of the population as a whole, in particular its size 
distri bution. This is akin to the c oncept of 'progenitor 
bias' ijvan Dokkum h Franxl [2"001h . As it turns out, a 
significant size evolution is expected even in the absence 
of size evolution of individual galaxies. 

Then we continue and include a prescription for 
subsequent merging based on the cosmologic al simu- 
lation s of the dark matter halo assembly by iLi et al.l 
()2007| ). The merger rates that we deri ve are in agree- 
ment with those f ound emp i rically (e.g.. lBell et al.l[2006l : 
iLotz et all l2008t iLin et al.l l2008t ). Because we exam- 
ine the merger history after the last, major episode of 
star-formation activity and the transformation into an 
early-type galaxy, b y definition, these mergers are dis- 
sipationless or 'dry' (jvan Dokkuml[200llBell et al.l[2006l : 
|Naab et al. 2006). Under reasonable assumptions for the 
relationship between the properties of progenitors and 
descendants based on numerical sim ulations of merging 
gas-poor progenitors (e.g., Ciotti fc van Albadal 
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^onzalez-Garci'a & van Albada 2003; Nipoti et al.l 12003 
Bpylan-Kolchin et al. 2005; Robertson et al. 200q), we 
can predict how individual galaxies grow. This evolu- 
tion of individual galaxies is then superimposed on the 
evolution of the population (Section |4]). A comparison 
with the observations allows us to investigate whether the 
continuous growth of the early-type galaxy population, 
with galaxies experiencing 'puffing' through dry mergers, 
is able to explain the observed size evolution with red- 
shift, or whether other, perhaps unknown, physical mech- 
anisms must play a significant role. Potential caveats, 
known problems, and testable predictions of our model 
are discussed in Sections [S] and El We adopt as cosmolog- 
ical parameters (r^M, f^A, cs) = (0.3, 0.7, 0.7, 0.9). 

2. SUMMARY OF HIGH-REDSHIFT OBSERVATIONS 

Before describing and testing our model, we summarize 
the relevant high-redshift measurements that we seek to 
explain. These are the evolutions of the characteristic 
size and the comoving mass density of early-type galaxies 
from z = 1 and z = 2 to the present. 

2.1. The Early-Type Galaxy Population at z ^ 1 

Size evolution between z ~ 1 and the pres e nt ha s 
recently been measured by Ivan der Wei et aD ((2008h . 



who found, at a given dynamical mass, i?off (z) oc (1 -|- 
^-)-o.98±o.ii fgj. galaxies with masses M > 10" Mq. 
This corresponds to SRcs{l) = i?off (l)/^cff (0.06) — 
0.54 ± 0.04 (z = 0.06 is the center of the redshift bin in 
which we select our sample of nearby early-type galaxies 
described in Section |3| and will be used throughout as a 
benchmark). 

The second observational quantity we compare our 
model with is the evolution of the comoving mass density 
p. We derive the change in the comoving mass density of 
early-type galaxies with masses larger than 10^^ Mq from 
the evolution of the normalization {(p*) of t he luminos- 
ity fun ction for red galaxies as measured bv lFaber et al.l 
()2007D . Even though not aU massive, red galaxies have 
early- type morphologies, the fraction of early- type galax- 
ies among th ese changes little between z = 1 and th e 
present (e.g., iBell et al.l[2004at Ivan der WeleTalllMl . 
We assume that the characteristic mass at the 'knee' of 
the mass function M* ~ 7.1 x 10^" Mq does not change, 
implying that the evolution in 0* is identical to that in p. 
This assumption should hold because the evolution of the 
M/L, as derived from evolution of the fundamental plane 
zero point and the observed evolution of L* cancel. How- 
ever, it cannot be ruled out at this stage that M* changes 
by ~ .1 dex. From the results presented bv lFaber et al.l 
(j2003), we derive that Sp{l) = p(l)/p(0.06) = 0.35±0.13 
for galaxies with mass M > 10^^ Mq, where the error 
includes an uncertainty of 0. 1 dex in the evolution of M* . 

2.2. The Early-Type Galaxy Population at z 2 

The observational constraints at higher redshifts are 
far less secure than at z < 1 due to limited sample 
sizes and systematic uncertainties. Simply extrapolating 
the z = 1 size evolution measurement mentioned above 
to z = 2, without regard to other measurements, gives 
(5i?off (2) = 0.36 ± 0.04. Compihng all available data sets 
at hig her redshifts (up to z = 2.5), Ivan der Wei et al] 
((200l found marginally faster evolution of Rcs{z) cx 



(1 + ^)- 



-1.20±0.12 



i.e., (5i?cff(2) = 0.29 ± 0.04. Assum- 



ing that the difference between the robust, but gener- 
ously extrapolated, z ~ 1 results, and the direct, but 
systematically uncertain, z ~ 2 results is indicative 
of the true, systematic error, we conservatively adopt 
(5i?ofr(2) = 0.3 ± 0.1 as the most realistic estimate of the 
size evolution between z = 2 and the present. We note 
that all measurements have been obtained for galaxies 
with masses M > 10^^ Mq. At lower masses samples 
are non-existent or severely biased. 

The observational constraints on (f)* and hence p 
are less secure at z 2 than at z ~ 1, again 
because of systematic unce rt aintie s. The latest es- 
timate by iMarchesini et al.l ()2008D yields 5p(2.4) = 
0.14lg'Qg for all galaxies with stellar masses larger than 
10^^ Mc.^^. Other recent deter minations ([Eisner et al.l 
l2008tlPerez-Gonzalez et al.ll2008D are consistent with this 



Kriek et al.l (|2008h have shown that already at 



estimate. 

z ~ 2.4 the majority (56%^^2%) such massive galaxies 
have red colors, and that those show no or little evi- 



The mass estimates in that study assume the lSaloetej p95l l 
initial mass function for which, given the other uncertainties, the 
difference between stellar mass and dy namical mass is nogligi- 
ble for the mass range of interest here l|van der Wei et al.i i2006l : 
IBorch et aLll2006i : ICappellari et al.ll200^ . 
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Fig. 1. — (a): The age-size relation for present-day early-type galaxies in nine different dynamical mass bins. The age is the light- weighted 
stellar population age. The distance between the same-colored error bars is equal to the scatter (Icr) in size. The median size and age in 
each bin, indicated by the stars, both increase with mass (from blue to red), reflecting the well-known mass-age and mass-size relations. 
However, at fixed mass small galaxies are older than large galaxies, (b): The age-size relation as in panel (a), but here in nine bins of 
stellar velocity dispersion instead of M^y^. The same overall trend is seen in the sense that galaxies with high cr* are, on average, larger 
and older than galaxies with low cr*. However, at fixed cr* the age-size trend differs importantly from that seen at fixed Mdyn^ age does 
not depend on size (or size does not depend on age) at fixed cr*. (c) and (d): The same as panels (a) and (b), respectively, but with H(5 
line strength instead of age along the y-axis. The trends are similar as those seen with age, albeit with larger uncertainties. Since H<5 is a 
robust age indicator, the similar behavior in this figure of H(5 and age demonstrates that the trends seen with age are not caused by model 
uncertainties in the age estimates, and cannot be artificial, (e) and (f): The same as panels (a) and (b), respectively, but with metallicity 
instead of age along the y-axis. As with age, there is no correlation between size and metallicity at fixed cr, . We note that the size of the 
SDSS spectroscopic aperture cannot explain the observed trends (see Section |3. II for details). 

of 17,483 nearby early-type galaxie s at redshifts 0.04 < 
z < 0.08 has been constructed by [Graves et al.l ()2008D 
from t he Sloan Digital Sky Survey (SPSS) database 
(DR6. lAdelman-McCarthv et all [20081 ). Only galaxies 
with measured velocity dispersions, on the red sequence, 
either without emission lines or with high [OII]-to-Ha 
ratios and with concentration parameters C > 2.5 are 
included in the sample. These criteria effectively exclude 
star-forming galaxies, but include genuine early-type 
galaxies with nuclear activity (see Yan et al. 2006). For 
the description of the velocity dispersions and the deter- 
mination of effective radii, as well as o ur dynamical mass 
estima tor (Mdyn oc i?fT^), we refer to Ivan der Wei et al.l 
(|2008D . In addition, we require that luminosity-weighted 
stellar ages are available from Gallazzi ct al. (2005); 
16,279 out of 17,483 galaxies have such an age estimate. 

3.1. Correlations between Global Properties of 
Early-Type Galaxies 

It is well known that early-type galaxies with large 
masses and velocity dispersions have older stellar popula- 
tions than those wit h low masses and velo c ity d ispersions 
(iTraeer et all 120001: iThomas et al.l 120051: [Gallazzi et all 



dence for star-formation activity. This number is some- 
what lower than for the nearby massive galaxy popula- 
tion which has a red sequence fraction of 77%. In addi- 
tion, despite the quiescent nature of the red z ^ 2 galax- 
ies, it is as yet unclear what their morphologies are. They 
are barely spatially resolved, such that it is impossible to 
tell whether they have smooth surface brightness profiles, 
and constraints on the shape of their surface brightness 
profiles are rather poor, even though iToft et all (|2007f ) 
showed that exponential profiles are perhaps a better 
representation than De Vaucouleurs profiles. Given these 
constraints on the evolution of the red/early- type frac- 
tion among the massive galaxies at z ~ 2, we arrive at 
Sp{2A) = O.lOlQog as our best estimate of the comoving 
mass density of early-type galaxies a,i z = 2.4 with mass 
M > 10" Mq. 

3. A SIMPLE EMPIRICAL MODEL I: 
DENSITY REFLECTS FORMATION EPOCH 

We develop an empirical estimate of the formation 
epoch of early-type galaxies based on the correlations 
between the global properties of the present-day pop- 
ulation of early-type galaxies. In this approach, indi- 
vidual early-type galaxies are assumed not to evolve af- 
ter their formation in either mass or size. A sample 



l2006riGraves et al.l |2007D. In Figure [fji we also show 
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Fig. 2. — Size vs. stellar population age for present-day early- 
type galaxies in a narrow range of velocity dispersion. The line with 
the errorbars is a running median, while the other lines indicate 
the scatter (Icr and la). The scatter is of order 10% and the 
residual from the running median does not correlate with any other 
parameter. 

that in our SDSS sample massive galaxies have older stel- 
lar populations. In the same figure, it can also be seen 
that, naturally, size increases with mass. As a corollary, 
we have that, loosely speaking, galaxies with large sizes 
have older stellar populations than galaxies with small 
sizes. 

For these broad, global trends it does not matter what 
tracer of mass is adopted. Dynamical mass and velocity 
dispersion both show similar correlations with size and 
age (compare Figures [1^ and[T]3). However, the differ- 
ence between the dynamical mass and the velocity dis- 
persion comes to light when we dissect the sample, and 
look at the relation between age and size at fixed dy- 
namical mass or at fixed velocity dispersion. At a fixed 
dynamical mass, there is an anti-correlation between size 
and age: large galaxies have younger stellar populations 
than small galaxies (Figure [T^). This trend persists over 
a mass range of almost two orders of magnitude, from 
~ 10^" Mq to ~ 10^^ Mq, and there is no clear indi- 
cation that the slope of the age-size relation depends on 
mass. 

The correlation between age and size at fixed mass 
implies that the zero point and scatter of the mass-size 
relation will evolve with redshift, even if we assume that 
individual galaxies do not change in size. This is simply 
due to the fact that larger, younger galaxies did not ex- 
ist yet at earlier ep ochs. This is, in essenc e, the model 
prediction made bv lKhochfar fc Silkl (|2006D . The result 
shown in Figure [1^ provide the first empirical evidence 
for this process. 

It is perhaps surprising that the picture changes quite 
dramatically if we look for an age-size relation for galax- 
ies with a given velocity dispersion. In Figure [T]d, it can 
be seen that at a given velocity dispersion there is no ev- 
idence that small galaxies are older than large galaxies. 
On the other hand, if one realizes that Mdyn oc -RcffC*, 
the size-age relation has to be flatter at fixed ti* than 
at fixed Mdyn (or may even be reversed). The fact that 
the relation should disappear altogether at fixed cr* is, 
however, not obvious. In Figure [21 we show for one par- 
ticular narrow bin around cr* = 200 kms~^ the size-age 
distribution. Again, there is no sign that age depends 
on size (or vice versa). Moreover, the scatter is remark- 



ably uniform and rather small (~ 10%). We could not 
identify a parameter that correlates with the residual in 
the size-age rela tion at fixed g^,. T his includes environ- 
ment: using the lYang et al.l ()2007f ) SDSS group catalog, 
wc find that the scatter in age does not correlate with 
group membership, the mass of the group, the distance 
to the group cente r, or whether galaxies are c entrals or 
satellites (see also, Ivan den Bosch eralll2008b[ ). 

The age estimates are unavoidably model dependent, 
and may therefore be hampered by systematic uncertain- 
ties. In Figure [Tt and[T]i we replace stellar population 
age by HJ absorption line strength. The strength of this 
Balmer line is relatively insensitive to metallicity and is 
therefore a good tracer of age. Although the random un- 
certainties are substantially larger, the picture remains 
the same: at fixed dynamical mass the H(5 line is stronger 
for larger galaxies, implying a younger stellar population, 
whereas at fixed velocity dispersion there is no correla- 
tion between size and line strength. We conclude that 
the trends with stellar population age in Figure [T^ and 
[T)3 are robust and not artificial. 

In Figures [T^ and [if, we show the correlation between 
size and metallicity. As with age, for the population 
as a whole there is a positive correlation between size 
and metallicity, as expected from the combination of 
the mass-metallicity relation and mass-size relations. At 
fixed velocity dispersion (Figure[T]:), however, there is no 
correlation between size and metallicity, which is consis- 
tent with the idea that the velocity dispersion traces the 
potential well depth. By necessity, we also find that at 
fixed dynamical mass, the metallicity is smaller for large 
galaxies (see Figure [TJj). 

The SDSS spectroscopic aperture does not sample 
the entire galaxy; therefore, some aperture effects may 
be expected: age and metallicity gradients could in- 
troduce artificial trends in Figure [T] However, aper- 
ture effects would cause a positive correlation between 
age/metallicity and size, not the observed anticorrela- 
tion. Therefore, we conclude that aperture effects do not 
strongly affect our analysis. 

The emerging picture is consistent with previous work. 
It has been shown that at fixed cr ^, there is no color- 
magn itude relation (jBernardi et al.l l2005t [Graves et al.l 
[20081 ). which indicates that velocity dispersion, more so 
than luminosity, determines the properties of the stellar 
population, unless some sort of conspiracy renders trends 
invisible. However, this possibility is excluded by the 
lack of a trend in Figures [TJd, [TJi, and [if. In addition, 
[Chang et al.[ ([2006f) find that, at fixed velocity dispersion, 
there is almost no correlation between stellar mass and 
absorption line indices. 

Very recently, several authors have noted similar trends 
(or lack thereof) as in our Figure 1. [Graves et all ([2009[ ) 
show that, at fixed velocity dispersion, age and metal- 
licity do not correlate with galaxy size (se e our Figures 
[If. and [If). Moreover. [Shankar fc Bernardj ([2009[ ) derive 
the same trend that we show in Figure la: at fixed mass, 
old galaxies are smal ler than young g alaxies. 

Furthermore, [Bell fc de Jond ([2000D and 
[Kauffmann et all ([2004[ ) showed that for any type 
of galaxy, stellar mass surface density is a better 
predictor than stellar mass of the star-formation history 
and current star-formation activity, and surface density 
is closely related to velocity dispersion. This was made 
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explicit by iKaufFmann et al] (|2006D . who showed that 
there is a sharp transition in the fraction of quiescent 
versus star-forming galaxies at particular values of both 
surface density and velocity dispersion. 

These previous results and our result that age and size 
at a given cr, are uncorrelated (Figure [1)3) tell us some- 
thing about the significance of (and perhaps the) 
fundamental characteristic of a galaxy. We now pursue 
the idea that a* can be used as a predictor of the epoch 
when star formation was truncated and an early-type 
galaxy emerged. 

3.2. The Formation Epoch of Early- Type Galaxies 

As a starting point, we take the observation by 
iKauffmann et al] ()2006( ) that at the present epoch the 
transition between star-forming and quiescent galaxies 
takes place at cr* ~ 125km s~^. Truncation of star forma- 
tion may be associated with the morphological transfor- 
mation to an early-type galaxy, and so we may suppose 
that newly formed early-type galaxies in the present-day 
universe have a velocity dispersion of ctet = 125kms~^. 
This transition is not infinitely sharp; in reality there are 
late-type galaxies with higher cr* and early-type galaxies 
with lower cr, . We ignore this scatter in our model, but in 
Section 15.11 we will briefly comment on the consequences 
for our model predictions. 

Given such a velocity dispersion threshold for star- 
formation activity in the present-day universe, we may 
expect that such a threshold existed at earlier epochs as 
well, but not necessarily with the same value. The re- 
lation between age and velocity dispersion (Figure [TJa) 
implies that uet would have had to be higher in the 
past. Observational evidence for the exis tence of ctet a t 
high redshifts was recently found by Franx et al.l (|2008D , 
who used data extending to z ~ 3 to show that there ex- 
ists a surface mass density threshold, Eet, above which 
star-formation activity drops suddenly and significantly. 
Interestingly, this threshold was shown to be larger at 
higher redshifts. We parametrize the velocity dispersion 
threshold for the truncation of star formation and the 
transformation into an early- type galaxy as follows: 

crET(^) = 125kms-i (1 + z)". (1) 

Consequently, given the velocity dispersion of a galaxy 
we can estimate the redshift at which its star formation 
was truncated and it attained its early-type morphology: 

Although precise constraints on the amount of evolu- 
tion are still lacking due to s ystematic err o rs an d small 
sample sizes, the results from iFranx et al.l ()20Q8[ ) imply 
that a must lie in the range l/4<a<3/4. In the 
following, we will adopt a = 3/4 as the standard value, 
which turns out to provide the best results. We will dis- 
cuss the dependence of our results on the chosen value 
and the range allowed by the data in Section O We note 
that the redshift dependence of (Tet is consistent with 
the long-standing result from fundamental plane stud- 
ies at high redshift that high-mass galaxies have stel- 
lar populations with high format i on re dshifts (z > 2, 
see Ivan Dokkum fc van der Marell l2007l and references 



therein), and also with the observation that this forma- 
tion redshift depends on mass, low-m ass galaxies having 
youn ger stellar popul ations epochs (jvan der Wei et al] 
2005: Treu et al.ll2005[ ). In addition, this is in agreement 
with archaeological studies of ne arby early-ty pe galaxies 
(e.g.. lThomas et al. 2005; ,Gallazzi et al.l[2008i) . 

To quantify the model predictions for earlier epochs, 
we apply Equation [1] to our nearby galaxy sample, only 
retaining those with ct, > aET{z) as galaxies that were 
early types already at redshift z. In Figure [3] we show 
how the early-type galaxy population would build up at 
different present-day masses for a — 3/4. Half of the 
galaxies with present-day masses exceeding 10^^ Mq had 
become early types by z ~ 2, while half of those with 
masses exceeding 10^^ M© had become early types by 
z 1. This behavior is reminisc ent of the pred ictions 
by, e.g.. lDe Lucia et all (2006) and lNeistein et al.l ([2006), 
who also predict that massive galaxies formed their stars 
earlier than less massive galaxies. 

In addition to comoving mass density evolution, Equa- 
tion[2]also implies size evolution at fixed dynamical mass: 
at a given mass, large galaxies have smaller velocity dis- 
persions and therefore became early-type galaxies more 
recently. Figure 2] shows the predicted evolution in the 
size distribution from the present to z = 1 and z = 2. 
The behavior of the model is straightforward: at progres- 
sively higher redshifts the minimum velocity dispersion 
of a galaxy (the dotted lines in Figure 2]) is higher, as 
defined by Equation [1] The resulting evolution in both 
the number of galaxies and their size distribution can 
readily be seen. How does this model compare with the 
observations described in Section [2]? For z = 1 , we find 
Sp{l) — 0.61 and (5i?off(l) = 0.81, both for galaxies with 
mass M > 10^^ Mq, and where the latter is the aver- 
age offset from the local mass-size relation, computed in 
log-space. The actually observed evolution is faster (see 
Section 12.11 and Figure [5|) , significantly for the evolution 
in average size and marginally for the mass density evo- 
lution. For z = 2 the model predicts 5p{2) = 0.09 and 
(5i?eff(2) = 0.71. The comoving mass density in this case 
agrees well with the observed value. The observed size 
evolution, however, is much stronger than this model pre- 
dicts (see Section [22] and Figure [5]). 

Despite the quantitative disagreements, it is encour- 
aging that this model, which only prescribes a redshift- 
dependent cr-threshold when galaxies stop forming stars 
and become early types, implies trends in the right di- 
rection. These results suggest that the continuous 'top- 
down' emergence of new early-type galaxies can explain 
roughly half of the observed size evolution, with indi- 
vidual galaxies not changing after their initial transfor- 
mation. However, an indication that indiv idual galaxies 
must change as well was recently found bv lTruiillo et al.l 
(|2009( l: the local abundance of galaxies with the same 
properties as the compact, high-redshift objects is very 
low, implying that such galaxies did not survive in their 
original form. It is clear that size evolution of individual 
galaxies has to occur as well. 

4. A SIMPLE EMPIRICAL MODEL II: 
GROWTH THROUGH MERGERS 

We now augment the simplistic model from Section [3] 
by a prescription to include subsequent merger activity 
of galaxies after the epoch of transformation into early 
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Fig. 3. — Predicted evolution with redshift of the comoving 
number density of early-type galaxies with different masses based 
on the velocity dispersion threshold specified in Equation |2] The 
horizontal dotted line indicates a fraction of 50%. High-mass galax- 
ies emerged as early-type galaxies earlier than low-mass galaxies. 
Note that this evolutionary picture does not include merging, as 
described in Section |4] 

types. After constructing merger histories, we describe 
how mergers affect, in our model, the masses, sizes, and 
velocity dispersions of early-type galaxies. 

4.1. Merger Histories 

After they become early types, galaxies still undergo 
subsequent evolution in terms of merging and/or mi- 
nor episodes of star-formation activity, without chang- 
ing their overall morphological properties, apart from, 
perhaps, brief periods of time after accretion events. 
That this can lead to substantial size evolution has been 
shown by Na ab et all (p007). To include such evolution 
schematically into our model we implement merger his- 
torie s based on simu lations of dark matter halo assembly 
(l20nl 

Li et all (pool , we define n as the inverse 
of the fractional mass increase in a merger. In other 
words, if the mass of the main progenitor is M^, then the 
mass of t he me rger remnant is Mf = (1 -I- l/n)Mi. The 
iLi et al.l ()2007f ) simulations demonstrate that the num- 
ber of mergers does not depend on n, i.e., the probability 
distribution of n is flat. To realize merger histories we 
can therefore assume a series of accretion events with ran- 
domly chosen n. We verify this simplified approa ch by at- 
tempt ing to reproduce the merger histories that iLi et al.l 
()2007t l inferred, starting with seed halos that have 1% of 
their present-day mass. When we allow mergers in the 
range 1 < n < 100, we construct an ensemble of merger 
histories with the same numb er of mereers w ith n < 3, 
n < 4, and n < 6 as found bv lLi et "aTl ()2007f ): see their 
Figures 11 and 12. In addition, our scatter in t he num- 
ber of such mergers is also close to that found bv lLi et all 
(|2007( ). We conclude that it is sensible to approximate 
the merger histories of dark matter halos by assuming a 
series of merger events with randomly chosen n in the 



from ILi et a 
Following 



range 1 < n < 100. 

Since we are not interested in the merger history of a 
halo growing from 1% of its present-day mass, but rather 
in the merger history of a halo between z — 1 ov z = 2 
and the present, we need to know by how much halos 
grow over that period of time. Contrary to the distri- 
bution of n, the grow th with time de pends on halo mass 
(|van den Boschll2002t IWechsler et al.ll2002D . We focus on 
halo masses in the mass range Mhaio = 10^^ — IO^^Mq, 
which corresponds to the halo masses of gala xies with 
masse s in the range M ^ 10^^ — IO^^Mq. In the lLi et"aLl 
(|2007h simulations, halos in this mass range had as- 
sembled ^ 50% of their present-day mass by z = 1 
(M^=i = M-io^o^oo^ where 0.09 is the scatter), and 



30% by 2 = 2 (M^=2 = M:^ 



-0.55±0.28 



Combining the elements discussed so far, we have 
the following practical method to construct the possi- 
ble merger history between z — 1 and the present for a 
halo with present-day mass M . We assume a series of 
mergers with mass ratios l/n^, with 1 < < 100 a se- 
quence of randomly chosen real numbers. For example, 
the most recent merger involved a main progenitor with 
mass riiM / [ni + 1) and a minor (or, accreted) progenitor 
with mass M/{ni + 1). The main progenitor is, in turn, 
assumed to be the product of a merger with ratio l/ri2. 
The main progenitor is followed in this manner, assum- 
ing a sequence of mergers with ratio l/rii, up until the 
point that the expected mass for the main progenitor at 
z = 1 is reached. This expected mass is Mz=i = M^^q, 
where x is drawn from a Gaussian distribution with a 
mean —0.30 and a standard deviation 0.09. To construct 
the merger history between z = 2 and the present, the 
only difference is that we choose x from a Gaussian distri- 
bution with mean —0.55 and a standard deviation 0.28. 
In practice, in order to conserve mass, the mass ratio 
of the final event is rounded such that the M^^i is pre- 
cisely M|Lg. In case x > 0, which we formally allow, we 
assume that no mergers have occurred. 

This process descr i bes th e merger history of a halo. 
Ivan den Bosch et al.l (|2008af l found that the vast major- 
ity of the galaxies in the mass range of interest here 
should be central galaxies. We, therefore, assume that 
the merger history of the halo directly corresponds to 
the merger history of the galaxy. We incorporate one 
exception to this rule into the model: mergers involv- 
ing halos with mass Mhaio < 10^^ M© do not increase 
the mass of the galaxy, but only the mass of the halo. 
This takes into account that such sm all halos contain 
a much smaller fraction of star s (e.g., ICoorav &: ShethI 
l2002Hvan den Bosch et~ani2007l ). The growth in galaxy 
mass is therefore somewhat slower than the growth in 
halo mass. As it turns out, our model predictions are 
not strongly affected by this. In Section 15. 2i we discuss 
the effect of possible delays between mergers between ha- 
los and mergers between their occupying galaxies. 

For each galaxy in the nearby sample, we construct a 
Monte Carlo realization of its merger history as described 
above, which predicts by how much the mass of the main 
galaxy grows between z = 1 (z = 2) and the present, and 
how the accreted mass is distributed over a number of ad- 
ditional progenitors at z = 1 (z = 2). Hence, we have 
two sets of merger trees: one that describes the z = 1 
progenitors and another that describes the z = 2 progcn- 
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Fig. 4. — Predicted reds hift evolution of the mass-size relation. The left-hand column of panels shows the predicted evolution for the 
mod el described in Section 13.21 i.e., without merging. The model presented in the middle column includes merging as described in Section 
14.31 only considering the main progenitors of pre sent -day early-type galaxies. The model presented in the right-hand column includes, in 
addition, all progenitors, as described in Section 14.41 The top row (in which all panels are identical) shows the observed, local mass-size 
relation, which is the starting point of our model, the gray-scale corresponding to the number of galaxies. The solid line is a running median, 
and the dashed lines indicate the scatter (Itr). These lines are repeated in every panel to guide the eye. The diagonal dotted line corresponds 
to (7, = 125 kms^^, which is the present-day velocity dispersion threshold, (tet> for the formation of early-type galaxies according to 
Equation[T] The middle row shows the predicted mass-size relation for z = 1. The dotted line now corresponds to ct* = 210 kms~^, which 
is o"Ex at 2 = 1. The lighter gray scale reflects the evolution in the number of galaxies. In addition, their size distribution is different. The 
bottom row shows the predicted mass-size evolution for 2 = 2, where (Tet = 285 kms~^. In our model, the effect of dry mergers is that 
galaxies move parallel to the dotted lines, i.e., a remains constant (see Section|4)l. 



itors. The average number of mergers per galaxy with 
mass ratios exceeding 1:2 is ~ 0.25 between z — 1 and the 
present (~ 0.55 between z = 2 and the present); the av- 
erage number of mergers with mass ratios exceeding 1:4 
is ^ 0.40 0.90). Note that whereas the dry merger 
rate per unit time per early-type galaxy was higher in 
the past, the dry merger rate per unit time and per unit 
comoving volume was smaller due to the strong evolu- 
tion in the comoving density of early-type galaxies. The 
merger rates in our model are consistent with or lower 



than the available observational constraints on the dry 
merger rate between z ^ 1 and the p resent (|Bell et alJ 
l2(M iLotTeTall [20081 : iLin et al.l[2008t ). 

4.2. The Effect of Dry Mergers on Galaxy Sizes 

A crucial part of the model is the effect of merg- 
ing on the sizes of galaxies. Because we are inter- 
ested in the merger history of galaxies that are al- 
ready early types, we focus on dissipationless (dry) 
merging. According to most numerical simulations the 
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remnants of gas-poor mergers of equal-mass progeni- 
tors situated in dark halos have roughly twice the size 
of thos e progenit ors (e.g.. ICiotti &: van Albadal 
Gonzal ez-Garcia fc van Albadall2003l: Nipoti et al.l 



2001 



2003 



Robertson et al.li2006l ): vet iBovlan- Kolchin et air(|2005l ) 



find that the size of the remnant is typically only 1.5 
times that of the progenitors. This difference may, in 
some cases, be explained by the assumed orbits: the 
bound orbit used by Boylan-Kolchin et al. may be re- 
sponsible for remnants that are smaller than found by, 
e.g., Ciotti et al., who use a parabolic orbit. However, 
not all simu lations with bound orb its lead to small rem- 
nants (e.g., [Robertson et al.l[2006f ). such that the rea- 
son for the discrepancy is not completely clear. In 
agreement with the majority of the predictions, we as- 
sume that equal-mass, dry merger products have dou- 
ble the size of the progenitors (see Section [S] for fur- 
ther discussion). The virial theorem then implies that 
the velocity dispersion remains constant, supported by 
the numerical simulations. Moreov er, t he results of 
iGonzalez-Garcia fc van Albadal (|2003f ) and iNipoti et al.l 
( 20031 ) suggest that for a dry merger with any mass ratio 
the size of the main progenitor increases linearly with 
mass, keeping its velocity dispersion constant. The im- 
plied decrease in surface brightness is such that merg- 
ers of this kind move galaxies roughly, but not precisely, 
along the fundamental plane (see, e.g. [Robertson et al.l 
[2006h . This is reassuring because scenarios in which this 
is not the case are difficult to reconcile with the observe d 
tightness of this scaling relation (e.g., iRix et al] Il994l ). 
These considerations allow us to directly and linearly link 
the growth in mass and size in our merger scenario. 

4.3. The Evolution of the Main Progenitors 

First, we investigate the evolution of the main progen- 
itors in the merger scenario described above, ignoring 
the smaller progenitors. Because we assume that merg- 
ers do not alter the velocity dispersion of the main pro- 
genitor, it is straightforward to combine the predicted 
merger history from Section 14.11 with the formation red- 
shift criterion described in Section 13.21 Figure |4] shows 
how the evolution in the number of early-type galaxies 
and their size distribution is altered when we augment 
the (Tet model from Section [3] by dry merging, only con- 
sidering the main progenitors. Relative to the scenario 
without merging, the number of high-mass galaxies de- 
creases more rapidly with redshift, while the number 
of lower-mass galaxies decreases less rapidly. The aver- 
age difference in size is predicted to be somewhat larger 
than in the absence of merging. The agreement with 
the observations improves considerably (see Figure [5]) : 
The merger model now predicts a redshift evolution of 
^i?off(l) = 0.68 and 6Res{2) = 0.45 for the average size 
at a given mass, and Sp{l) = 0.38 and Sp{2) = 0.03 for 
the comoving mass density. The agreement with the ob- 
servations is not perfect though; the observed evolution 
of the mean size is still ^ 30% faster than predicted. 

4.4. The Properties of Accreted Progenitors 

At redshifts z < 1, when few major mergers occur and 
the growth of halos is relatively slow, it is probably suffi- 
cient to consider only the main progenitors, as we did in 
the previous section. However, at higher redshifts, when 
major mergers are more frequent and halos grow fast, this 



simplification may become inappropriate. It is, therefore, 
also necessary to include the smaller progenitors in the 
analysis and assess the question how their masses and 
sizes affect the evolution of the population averages. In 
the following we do not consider merger activity of the 
smaller progenitors before their accretion onto a larger 
object; thus, some of the uncertainty remains. 

We make two assumptions regarding the properties of 
the minor progenitors. First, because we are mainly in- 
terested in dry merging, we assume that all progenitors 
are early-type galaxies. Second, we assume that the ve- 
locity dispersions of all progenitors are equal to the dis- 
persion of the main progenitor, and therefore also equal 
to that of the final descendant (see Section 14. 2p . The 
second assumption is seemingly quite extreme and is ob- 
viously not valid for minor mergers. However, since the 
goal is to predict only the abundance and sizes of pro- 
genitors with mass M > 10^^ Mq, this is not necessar- 
ily a problem. Moreover, Equation [1] implies that, in 
our model, all early-type galaxies at any epoch have a 
minimum velocity dispersion, ctet- At redshifts z ^ 2, 
where smaller progenitors matter, the number of early- 
type galaxies is a rapidly declining function of veloc- 
ity dispersion such that most galaxies available for dry 
merging will have comparable velocity dispersions. Thus, 
given the characteristics of our model so far, it is defen- 
sible that galaxies that partake in mergers at z ~ 2 with 
mass ratios 1:3 or 1:4 have similar velocity dispersions. 

The population of smaller progenitors by definition 
have smaller masses than the population of main pro- 
genitors. Their inclusion will mostly affect the predicted 
comoving mass density evolution at z = 2 (see Figure 
111). Quantitatively, the model predicts 6p{l) = 0.45 and 
Sp{2) = 0.05, both consistent with the data (see Fig- 
ure |S|). The effect on the predicted evolution in average 
size at fixed mass is less pronounced: the model predicts 
(5i?off(l) = 0.63 and 5Rcs{2) = 0.40, not much slower 
than observed (the remaining disagreement is of order 
1 - 2a). 

Overall, it is quite remarkable that such a simplistic 
estimate of the redshift at which galaxies transform into 
early types and the simple assumptions regarding their 
subsequent merger activity can come so close to repro- 
ducing the observations. 

5. MODEL UNCERTAINTIES 

The model presented in Sect ions |3| and |4| is a drastically 
simplified description of the formation and evolution of 
early-type galaxies. In this section we describe some of 
the uncertainties that arise from the simplifications and 
their potential impact on the model predictions. 

5.1. Uncertainties in the Estimate of the Formation 

Epoch 



As stated in Section 13.21 the value of the exponent 
in Equations [T] and m was chosen a posteriori, a = 3/4 
producing the best results. The model predictions are 
quite sensitive to the choice of a. Trying a range of values 
for a, we find that, for the 'standard' version of our model 
as described in Section 31 in order to be consistent with 
all data at the 2cr level, the exponent has to have a value 
in the range 0.65 < a < .85. T his is consistent with 
the results of iFranx et al.l (j2008f ). For lower values of 
a we find, in particular, that the predicted evolution of 
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Fig. 5. — Comparison of the observations described in Section 
[2] (points with error bars) with the model predictions for redshift 
evolution of the average size and fixed mass (the top panel) and 
the total, comoving mass density (the bottom panel) for all early- 
type galaxies with mass M(z) > 10" Mq. The dotted, dashed, 
and solid lines refer to different versions of our model, as in Figure 
[4] T he dotted line represents the model without merging (Section 
13.21 1: the dashed line represents the mo del with merging, only con- 
sidering the main progenitors (Section l4.31 : and the solid line rep- 
resents the standard version of our mod el, including merging and 
considering all progenitors (Section 14.41 1. Our dry- merging model 
successfully reproduces both size and mass density evolution from 
z = 2 to the present. 

the mean size is slower than observed. For higher values 
the predicted comoving mass density at z ~ 2 is lower 
than observed. Clearly, if we change the prescription for 
merging (see Sections 15.21 and | 5 . 3p . the allowed range for 
a also changes; however, in order for the model to remain 
consistent with the observed evolution of the comoving 
mass density, a has to be close to 0.75. 

A related problem is a possible systematic uncertainty 
in the velocity dispersions of nearby early-type galax- 
ies. This remains to be an issue , especially for the most 
massive galaxies ()Bernardill2007() . Even a 5% systematic 
error has major consequences for the z = 2 predictions 
because the number of galaxies with high dispersions de- 
clines very rapidly with increasing dispersion. This im- 
plies that the formation-redshift distribution is rapidly, 
and perhaps artificially, truncated as well. Increasing all 
velocity dispersions by 5% leads to a higher comoving 
mass density at 2; = 2 by as much as a factor of 3 (while 
the increase at z = 1 is just 10%), whereas size evolution 
is only slowed down by less than 5% over the entire red- 
shift range. We conclude that the predictive power of our 
model regarding the comoving mass density evolution at 
z ~ 2 is likely limited by systematic uncertainties in the 
velocity dispersions of massive nearby galaxies. 

A further simplification in our model is that there is 
no scatter in the relationship between velocity disper- 
sion and formation redshift. As a result, the predicted 
comoving mass density evolution will be overestimated, 
particularly at z ^ 2 where the model predictions depend 
strongly on the steepness of the mass function. It is be- 
yond the scope of this paper to properly implement such 
scatter, as this would require a more advanced synthesis 



between the formation epoch and merger history of the 
galaxies in our simulated samples. 

Finally, we assume that the truncation of star for- 
mation coincides with a morphological transformation. 
Even though those events may have the same cause, they 
need not occur simultaneously. 

5.2. Uncertainties in the Estimated Merger Histories 

We make the approximation (see Section HTTj) that the 
merging of dark matter halos is followed by the merging 
of the galaxies within Ai << <Hubbic due to dynamical 
friction. In reality, galaxies occupying merging halos may 
merge not at all, such as in the case of infall of a galaxy- 
sized halo into a cluster-sized halo, or they may merge 
quickly, such as in the case of two galaxy-sized, equal- 
mass halos. Generally speaking, the larger the mass- 
ratio of a halo-halo merger, the l onger the delay of the 
galaxy-galaxy merger (see, e.g ., iBovlan-Kolchin et al.l 
I2008t iKitzbichler fc Whiti 12008?) . It may, therefore, be 
necessary to invoke a delay for galaxy mergers with re- 
spect to halo mergers. To obtain an upper limit, we test 
the effect of a long, 2.5 Gyr delay on our model predic- 
tions, which is of the same order as the Hubble time at 
high redshift, and amounts to a few times the dynamical 
time scale of large clusters. In practice, this means that 
the galaxy merger activity between z — \ and the present 
reflects halo merger activity between z = 2 and z — 0.2. 
As a consequence, major merger activity (with mass ra- 
tios of 1:2 or less) is boosted substantially from 0.25 to 0.5 
mergers per galaxy between z = 1 and the present. This 
results in the stronger size evolution ((5i?cff(l) = 0.57 
and (5i?cff(2) = 0.35) which would be in excellent agree- 
ment with the data (see Figure [5]) . Both the observa- 
tional and model uncertainties are too large, however, to 
draw the conclusion that a delay in galaxy merger activ- 
ity is preferable, especially because the mass dependence 
of this phenomenon is not taken into account. 

A merger delay for the galaxies under consideration 
may indeed not be of great importance, as cosmologi- 
cal simulations that explicitly treat galaxy mergers as 
separate from halo mergers find si milar merger rates 
for galaxies as we fi nd for halos (jMaller et al.l l2006t 
iKhochfar fc Siigi20?)8l ). This, along with our finding that 
any merger delay only increases the redshift evolution of 
the mean sizes (i.e., our model in this sense provides a 
lower limit), implies that our conclusions are not compro- 
mised by the simple assumption that galaxies mergers are 
instantaneous refiections of halo mergers. 

5.3. Uncertainties in the Properties of Progenitors 

Our modeling assumes that in a dry merger, the size of 
the main progenitor grows proportionally to its mass in- 
crease. We justify this choice in Section^ but it is useful 
to discuss its effect on the model predictions. If a sub- 
stantial amount of orbital energy were transferred from 
the stellar component to the dark halo, the descendant 
would be more tightly bound than the progenitors. This 
would result in a fr actional size change les s than that in 
mass, which is what iBovlan-Kolchm et aTl (l2005l) find in 
t heir simulations. The results from lBovlan-Kolchin et all 
(200^) in fact suggest that dry mergers move galaxies 
along the slope of the present-day mass-size relation, in 
which case no redshift evolution in size at fixed mass due 
to merging would be expected. Our approach implies 
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that no energy is transferred to the halo, resulting in 
large, relatively loosely bound merger remnants (with an 
unchanging velocity dispersion) , in agreement with most 
other simulations (see Section [4. 2[) . This is why mergers 
play an important role in shaping the early-type galaxy 
population in our modeling. Conversely, we would infer 
that the strong redshift evolution of the mean size im- 
plies that the amount of energy transferred to the halo 
must be small. 

The above arguments make it very clear that our model 
would not work without the assumption that galaxies, 
at least the main progenitors, grow in size more or less 
proportional to their mass. The properties of the smaller 
progenitors before merging with the main progenitors are 
a separate, but related issue. In Section we justified 
our assumption that all progenitors partaking in a merger 
have nearly the same velocity dispersion. This breaks 
down at lower redshifts, where a reservoir of galaxies 
with a large range in mass (and dispersion) is available 
for merging. It is therefore useful to test the sensitiv- 
ity of the model predictions on this assumption. To do 
so, we consider the following alternative scenario: the 
smaller progenitors follow the present-day mass-size rela- 
tion with respect to the main progenitor, i.e., we assume 
that t heir size ratios are (Afaccr/Afmain)"'^^ (|Shen et al.l 
120031 : Ivan d cr Wcl ct al. 2008). The virial theorem re- 
quires lower velocity dispersions for the smaller progen- 
itors than for the main progenitor. By comparing this 
lower velocity dispersion with the threshold associated 
with the particular redshift that is being simulated (i.e, 
z = 1 or z = 2), we can decide in a very natural man- 
ner within the context of our model whether the smaller 
progenitor is a late-type galaxy or an early-type galaxy. 
If it has such a low dispersion that the implied formation 
redshift according to Equation [5] is lower than the simu- 
lated redshift, it is omitted from the sample, otherwise it 
is retained. Despite the substantially different approach, 
this exercise has a limited effect on the model predic- 
tions. The predicted size evolution is ~ 10% slower, and 
the predicted comoving mass density evolution ~ 10% 
faster. We conclude that the assumptions regarding the 
properties of the smaller progenitors are not of crucial 
importance to the model predictions. 

It is, in principle, possible that in a (minor) merger 
the relative increase in size is larger than the relative in- 
crease in mass. In the case of virialized progenitors and 
remnants, homology, and zero energy transfer between 
the stars an d the dark matter halo, follo wing the analyt- 
ical work of lBovlan-Kolchin et al.l ()2005[ ). we derive that 
the ratio of the radii of the merger remnant (Rr) and the 
main progenitor (i?m) is 

Rr_ ^ + 

where n, as before, is the mass ratio of the main pro- 
genitor and the accreted progenitor, and a is the slope 
of the mass-size relation for the progenitors {R oc M"). 
For n oo, we have that Rr/Rm — > 1 + 2/n (for a < 1), 
implying that for minor accretion events, the relative 
increase i n size is close to twice the increase in mass 
see also. iNaab eFal] [20091) . This may provide an ex- 
planation for the evolution of the most compact galaxies 
in the observed z > 2 samples (e.g., IZirm eTall I2OO7I: 



Ivan Dokkum et all 12008) . which our model does not re- 
produce. The difficulty in the context of our model, 
however, is that knowledge of the nature of very small 
accreted galaxies, which in this scenario would be the 
main contributors to size evolution, becomes essential. 
It seems unrealistic to suppose that low-mass galaxies 
are not gas rich and are not forming stars. As a con- 
sequence, mergers will be dissipative, with a different 
effect on the evolution of the scaling relations (see, e.g., 
'Robe rtson et al.ll2006t ICiotti et al.ll2007t iHopkins et all 
i2009bl ). A more complete description of the merger his- 
tory of galaxies, including a prescription for dissipation, 
will be required to fully explore the role of mergers with 
large mass ratios. However, if energy transfer to the halo 
is a minor factor, the large number of expected minor ac- 
cretion events could be a viable explanation for the most 
compact, observed galaxies at 2 > 2. 

If minor mergers that are gas poor dominate the 
growth of early-type galaxies, size evolution can be 
regarded as the assembly of a low-density envelope 
around a high-density central region. In such a sce- 
nario, the compact, high-redshift galaxies survive as 
the centers of present-day, massive galaxies. The prop- 
erties of the centers of local early-type galaxies and 
the distant compact galaxi es are, i ndeed, not signif- 
icantly different (Bezanso n et al.l 12009: Hopkins ct aQ 
l2009aD . On the other hand, there is no indication that 
the outer envelopes of early-type galaxies, the wings 
of De Vaucouleurs-ty pe surface brightness pr ofiles, have 
evolved since z 1 (Ivan der Wei et al]l2008h . We con- 
clude that it is thus far unclear, from an empirical per- 
spective, what the relative contributions of major and 
minor merging are. 

Summarizing, there are two elements of the model 
that are essential to any success in explaining the ob- 
servations: first, the evolving velocity dispersion thresh- 
old which is associated with the formation of early-type 
galaxies; second, the growth of the main progenitor of 
present-day early-type galaxies by means of the accre- 
tion of mass through dry mergers in such a manner that 
size grows proportionally to mass. 

6. TESTABLE PREDICTIONS 

Our model predicts mass-dependent size evolution, 
with the high-mass end of the mass-size relation chang- 
ing least with redshift (see Figure . This is mainly 
due to the higher formation redshifts of more massive 
galaxies. Such a trend is furthermore strengthened by 
the presumed properties of the accreted progenitors de- 
scribed in Sect ion |4^ This prediction is exactly opposite 
to the prediction by Khochfar & Silk (2006) that massive 
galaxies will display the strongest evolution. This seem- 
ing discrepancy may be explained if the s t ar-for mation 
histories of galaxies in the lKhochfar &: Silkl (120061 ) model 
are closely linked to the assembly histories of their halos. 
We explicitly separate these processes by invoking a cri- 
terion for galaxies transforming into early types that is 
unrelated to the assembly history, only considering fur- 
ther assembly through dry merger events after the epoch 
of their initial emergence as early-type galaxies. 

Unfortunately, the current observations are not of 
sufficient quality to decide this matter. The z 
1 mass-radius relation of the sample of early-type 
galaxies with dynamical mass measurements used by 
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Ivan der Wei et all ()2008[ ) has a somewha t steeper slope 
than the local mass-size relation (see also iFerreras et alJ 
|2009[) . but this is a marginal effect (of order 2a), and 
which, in addition, could also be explained as a differ- 
ence between cluster and field galaxies. Interestingly, 
the z ~ 2 samples show a hint th at the reverse may 
be the case ()van der Wei et aLll2008f ). with samples with 
the largest masses showing the largest offsets from the 
local mass-size relation, in qualitative ( but not quantita- 
tive) a greement with the predictions bv lKhochfar fc SiI3 
(gOOi). 

Another feature predicted by our model is a decrease 
with redshift in the scatter in size at fixed mass (see 
Figure [3]) . This is simply because at fixed mass the 
larger galaxies are younger: at earlier epochs the larger 
galaxies at a given mass were not yet members of the 
early-type galaxy population, such that the scatter was 
smaller than it is nowadays. How rapidly the scatter is 
predicted to evolve with redshift depends on the imple- 
mentation of merger activity, but regardless of the de- 
tails a decrease with redshift is expected. Therefore, a 
measurement of the scatter in the high-redshift mass-size 
relation mainly tests our idea that there is a relation be- 
tween age and size at fixed mass, which we already know 
to exist given the archaeological properties of present-day 
early-type galaxies (see Figures [2] and [1]) . The effect of 
dry mergers on the scatter is, in our model, of secondary 
importance. The z ~ 1 observations are consistent with 
a non-evolving scatter (within 2a), but the best-fitting 
value points toward smaller scatter at earlier epochs (see 
Ivan de r Wcl et al. 2008), in agreement with our model. 
Clearly, larger samples are required to confirm this ten- 
tative result. 

Our study here concerns several fundamental scaling 
relations for early-type galaxies. It is therefore interest- 
ing to consider the implications of our model for another 
fundamental scaling relation that we have not mentioned 
so far: the connection between s uper-massive black holes 
and t heir host galaxies (e.g.. iKormendv fc Richstond 



1S95'; 'Magorrian et al.' '1998': Ferrarese fc MerrittI 120001 : 
Gebhardt et al. 2000; Graham et al. 200lT). In our dry- 
merger scenario, the growth of galaxies is obviously con- 
nected to the growth of the mass of the central black 
hole (Mt). We assume that cr* remains constant in dry 
mergers (Section 14. 2p . whereas the mass of the halo, the 
galaxy, and that of the black hole grow more or less pro- 
portionally. Hence, the evolution of M, is more closely 
linked to the evolution of the total mass of both the halo 
and the galaxy than to a^. However, note that this is not 
an argument against a fundamental relationship between 
(T, and M, that may have determined the initial shape of 
the scaling relations between black-hole and bulge prop- 
erties. It is as yet unclear what the intrinsic scatter in 
these scaling relations is, or whether one of the relations 
is fundam ental in the sense that the intrinsic scatter i s 
zero fe.g.- lTremaine et al.|[2002l : iFerrarese fc Fordll2005D . 
Nor is it clear if and in which direction the M, scal- 
ing relations evolve with redshift for non-active galax- 
ies, desp ite recent progress for galaxies with ac tive nu- 
clei fe.g.. lSalviander et al.|[2007l: iWoo et al.ll2008D . In the 
context of our model, evolution with redshift is expected, 
and the intrinsic scatter is unlikely to be zero. Scatter 
and redshift evolution in the Af.-cr, relation is expected 
due to continuous dry merging and the variety of the 



merger history of galaxies with a given ct, (see Section 
14. 1|) . The scatter should then correlate with the size of 
the galaxy (bulge). Furthermore, assuming that ct* and 
M, are initially related, scatter and evolution in the M,- 
Mdyn (bulge mass) relation is expected due to the ct*- 
dependent formation redshift of galaxies, and the scatter 
in the relationship between cr* and Mdyn (see Section 
13. 2|) . The scatter should then correlate with size and the 
age of the stellar population. However, the current ob- 
servational constraints on both the scatter of the present- 
day black-hole scaling relations and their evolution with 
redshift are not sufficiently accurate to provide evidence 
against or in favor of our model. Indications that the 
scatter indeed co rrelates with a third p arameter have 
been found (e.g.. [Marconi fc Hunfj [20 03^) , but the sta- 
tistical significance of these results is q uite marginal and 
their interpretation still debated (e.g.. lGrahamll2008l ). 

In addition to the M,-cr* relation, the model presented 
here also predicts some characteristic s of other scaling 
relations, such as the Faber-Ja ckson (|Faber fc JacksonI 
I1976D and M52 - cr, (Dressie r^ al.l I1987D relations. 
Equal-mass, dry mergers have little effect on the Mg2 — 
cr* relation, as both the metallicity and the velocity dis- 
persion remain unchanged. The scatter may decrease 
somewhat, as merging will lead to regression to the 
mean metallicity at fixed velocity dispersion. Unequal- 
mass mergers will slightly decrease the metallicity of 
the merger remnant compared to the metallicity of the 
main progenitor (because of the mass-metallicity rela- 
tion), which may somewhat increase the scatter in the 
Mg2 — <J* relation. If we assume that luminosity is a 
proxy for mass, our model predicts scatter in the Faber- 
Jackson relation. This reflects the variety in the merger 
history of galaxies and is equivalent to the expected scat- 
ter in the M.-cr* relation. Hence, our model predicts a 
correlation between the scatter in the Af.-cr* relation and 
the Faber-Jackson. Both of these should also correlate 
with galaxy size. The correlation between the scatter in 
Faber-Jackson relation and galaxy size is, of course, well 
known. The existence of the fundamental plane is, there- 
fore, consistent with a scenario in which dissipationless 
merging plays an important role. 

7. SUMMARY AND CONCLUSIONS 

We present a simple, empirically motivated model that 
simultaneously predicts the evolution in the mean size 
and the comoving mass density of early-type galaxies be- 
tween z = 2 and the present. A large sample of nearby 
early-type galaxies extracted from the Sloan Digital Sky 
Survey serves as the starting point. The redshifts at 
which galaxies transform into early-type galaxies (i.e., 
their 'formation' redshifts) are estimated based on an 
evolving velocity dispersion threshold, uet, above which 
galaxies have low specific star-formation rates. This is 
motivated in part by the observation that the stellar 
population age of present-day early-type galaxies is a 
simple function of their velocity dispersion. In addi- 
tion to a prescription for estimating formation epochs, 
several merging scenarios between the formation epoch 
and the present are explored. Merger trees are generated 
such that they match the results from simulations of the 
mass assembly history of dark matter halos. Assuming 
dissipationless ('dry') merging and a simulation- inspired 
prescription for the effect of such mergers on the size 
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of galaxies, we successfully reproduce both the observed 
evolution in the mean size and in comoving mass den- 
sity of early-type galaxies with mass M > IQ^^ Mq be- 
tween z = 2 and the present. Our model quantitatively 
explains the recently measured, substantial size evolu- 
tion of early-type galaxies. Our model does so by com- 
bining the continuous emergence of early-type galaxies, 
starting with systems with the highest velocity disper- 
sion, with subsequent merging. Recently, and from the 
very different perspect ive of hydrodynamical simulations, 
ins et al.l ()2009bf ) also considered such a combina- 
tion of growth of the population and growth of individual 
galaxies, and arrived at similar conclusions. 

We emphasize that our model for estimating the 'for- 
mation' redshifts of early-type galaxies is entirely phe- 
nomenological. We have, so far, not addressed the ques- 
tion why the velocity dispersion would be closely related 
to the formation redshift. Surely, this is related to the 
formation of bulges, perhaps regulated by feedback mech- 
anisms. An obvious connection can be made between the 
crucial role that velocity dispersion plays in our model 
and the importance of black hole growth. It will be in- 
teresting to see whether a theoretical framework can be 
constructed that explains from first principles that early- 
type galaxies which form during a certain epoch have 
velocity dispersions that are independent of mass. 

The most troubling issue of our model, which affects 
the merger part, is the effect of mergers on the sizes 
of the descendants. As mentioned in Section 14.21 some 
confusion remains in the literature describing the results 
from numerical simulations. Even if this issue is resolved, 
it will remain unclear how a multitude of minor accretion 
events will affect size evolution. In principle, a stronger 
effect on the sizes of the main progenitors than we assume 
in our model is physically possible. This may explain the 
presence of very compact galaxies in the z > 2 samples; 
however, dissipation will have to be incorporated in our 
model to make realistic predictions, which is beyond the 
scope of this study. 

Besides the possibly underestimated effect of minor ac- 
cretion events in our model, there are other potential 
solutions to the puzzle that the most compact galaxies 
present. First, the systematic uncertainties in the ob- 
servations may still be underestimated; the ultimate test 
would be a direct measurement of the velocity disper- 
sio n. Initial steps in this dir ection were recently taken 
bv ICenarro fc Truiillol ()2009l ). They find, based on an 
analysis of stacked spectra, that the velocity dispersions 
of quiescent galaxies at z ~ 1.5 may be somewhat higher 



than those of present-day galaxies with the same stel- 
lar mass. Second, physical mechanisms besides merging 
may play a role. A recently proposed idea is that quasar 
feedback is perhaps responsible for the observed size evo- 
lution, through the ejection of gas, simultaneously pro- 
viding an explanation fo r the cessation o f star formation, 
and the increase in size ()Fan et al.ll2008[ ). An indication 
that this cannot be the mechanism that fully explains 
the observed size evolution is that a substantial frac- 
tion of this evolution took place between z = 1 and the 
present; over this period, massive early-type galaxies are 
known not to undergo phases of intensive nuclear activ- 
ity. On the other hand, quasar feedback could contribute 
at higher redshifts, providing an explanation for the re- 
maining discrepancy between our model predictions and 
the most compact galaxies found at z > 2. 
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